'Local' Quantum Criticality at the Nematic Quantum Phase Transition 
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We discuss the finite temperature properties of the fermion correlation function near the fixed 
point theory of the charge nematic quantum critical point (QCP) of a metallic Fermi system. We 
show that though the fixed point theory is above its upper critical dimension, the equal time fermion 
correlation function takes on a universal scaling form in the vicinity of the QCP. We find that in the 
quantum critical regime, this equal-time correlation function has an ultra local behavior in space, 
while the low-frequency behavior of the equal-position auto correlation function is that of a Fermi 
liquid up to subdominant terms. This behavior should also apply to other quantum phase transitions 
of metallic Fermi systems. 
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A large number of metallic strongly correlated systems 
exhibit the phenomenon of quantum criticality. Typi- 
cally the quantum critical regime is accessed either by 
continuously tuning an external control parameter (e.g. 
magnetic field or pressure), or by chemical means (e.g. 
doping). This phenomenon has been observed, among 
others, in heavy fermion systems, itinerant ferromagnets, 
ruthcnatcs, in CDW systems, in high T c cuprate super- 
conductors, and other complex oxides. 

The theory of quantum phase transitions in metallic 
systems goes back to Hertz's seminal workQ and its 
subsequent development by MillisQ- It yielded a clas- 
sification of these QCPs closely analogous to the theory 
of equilibrium classical critical phenomena. It is an or- 
der parameter theory which focuses on the behavior of its 
collective modes, with the added essential physical effects 
of gapless fermionic excitations in the quantum dynamics 
and as a mechanism for dissipation. A key assumption of 
this theory is that the fermions, although gapless, con- 
tribute primarily to dissipation mechanisms (e.g. Landau 
damping) and determine the effective dynamic critical ex- 
ponent z. Following this logic one expects that at a QCP 
the system should be scale invariant and the behavior of 
all its observables be controlled by some fixed point the- 
ory. Although this is indeed correct, the actual behavior 
of the fermionic quasiparticles (and their very existence) 
has remained a major unsolved problem. Naively one 
would have expected that the fermion correlators should 
obey scaling as well, perhaps with a non-trivial anoma- 
lous dimension, as in one-dimensional Luttinger liquids. 
We shall see that this is not the case for dimension d > 1. 

Experiments in quantum critical metallic systems show 
that these QCPs exhibit a phenomenon known as 'lo- 
cal' quantum criticality Q. There is strong experimen- 
tal evidence for scaling behavior in both frequency and 
momentum, with a characteristic dynamic scaling expo- 
nent z, of the correlation functions of the relevant order 
parameter and its associated susceptibilities, and of all 
other correlations in the bosonic sector of these systems. 
Thus, the spin correlation functions measured by neutron 
scattering and the current and density correlation func- 



tions measured by light scattering, do exhibit scaling. In 
contrast, at these same QCPs, the fermionic correlators 
(measured in angle resolved photoemission (ARPES)) 
typically exhibit scaling in frequency but show essentially 
no scaling in momentum, thus the term 'local quantum 
criticality'. The explanation of these seemingly disparate 
behaviors has been a long standing puzzle. 

In this paper we develop a fully solvable fixed point 
theory of the of the quantum phase transition between 
a Fermi liquid (FL) and a charge nematic Fermi fluid 
(N) 4]. A charge nematic Fermi fluid is a translation- 
ally invariant state of a metallic system which breaks 
spacial rotational invariance (or the point group sym- 
metry of the lattice) spontaneously^. A charge ne- 
matic state has been seen in experiments in quantum 
Hall systems 0,0 (a QCP has not yet been observed). 
There is also convincing evidence for a charge nematic 
phase in the ruthcnate Sr3Ru207 (including quantum 
phase transitions) ■ It has also been suggested that 
the 'hidden order' phase of the heavy fermion compound 
URu2Si2 may be due to a p-wave analog of the charge 
nematic state, ^lj We will show here that the fixed point 
theory of the nematic QCP exhibits the hallmark of lo- 
cal quantum criticality: anisotropic scaling (with z = 3) 
for the bosonic excitations, while the fermions show scal- 
ing in frequency and no scaling in momentum (up to 
the effects of irrelevant operators). In fact, at finite 
temperature, we find that the behavior of the equal- 
time fermion correlator is strikingly similar to the re- 
sults found in Ref. for vortex- like operators in the 
(seemingly unrelated) quantum Lifshitz model [Til]. The 
concept of local quantum criticality, has been the focus 
of recent work using perturbative methods|l4|. We will 
attack t his p roblem usin g n on-perturbative bosonization 
methods [lU HfiL IfR 0, El Both approaches agree in 
the perturbative regime. 

The (Pomeranchuk) transition to a charge nematic 
Fermi fluid is a quantum phase transition in which the 
shape of the Fermi surface (FS) acquires a spontaneous 
quadrupolar distortion, and breaks rotational invariance 
spontaneously down to 7r-rotations. It is the simplest 
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example of an electronic liquid crystal phase[5|- It has 
dynamic critical exponent z = 3 so that in two dimen- 
sions, the fixed point theory is above its upper critical di- 
mension. A simple computation to one-loop in the RPA 
interaction done in Ref. |J] obtained a fermion self-energy 
lmT,(q = qf,u) cx uj 2 ^ 3 >w at low energies. By a scal- 
ing analysis, the interactions are thus relevant by a power 
of 1/3. Hence, an understanding of the fermionic proper- 
ties of this QCP requires a non-perturbative treatment. 

We recently used 1 131 the method of high dimensional 
bosonization[li 011 E |U E3 to develop a non- 
perturbative fixed point theory of the nematic order pa- 
rameter and of the fermion propagator. We demon- 
strated explicitly that this is a non-Fermi liquid with 
a fermion residue Z which vanishes at the QCP, while 
the single particle density of states, N(ui), remains es- 
sentially unchanged from its FL value. We also showed 
that the bosonization result for the fermion propagator 
reproduces the perturbative results of Refs.JJ, |22j, and 
showed that it treats all rainbow diagrams and vertex 
corrections on an equal footing. 

In this letter, we discuss its finite temperature behav- 
ior. We show that in the quantum critical regime and at 
the level of the fixed point theory, the equal-time fermion 
correlations are ultra local in space due to an infrared di- 
vergence which wipes out all correlations in space to all 
distances. This divergence is controlled by dangerous ir- 
relevant operators which then set the correlation length 
of the fermions. In the nematic phase, symmetry pro- 
tects the Goldstone modes at finite temperatures leading 
to the conclusion that the spacial correlation length of the 
fermions vanishes in the thermodynamic limit. Despite 
this extreme behavior, the correlations at equal-positions 
and at long times remain well behaved and the single par- 
ticle density of states, N(u>), is essentially unaffected. We 
also show that the equal-times function obeys a scaling 
form and utilize this behavior to describe the collapse of 
the FL as a critical phenomenon. 

We take as our starting point, a low energy effective 
Landau model of two-dimensional spinlcss fermions with 
a quadrupole density interaction: 



n = 



|q|-fc F <A 



(2tt)s 



^+f 2 (q)Q(q)-Q(-q) , (1) 



iV(0)/2(0) = F 2 < —1. The behavior that we discuss here 
should also apply to lattice systems at their QCP |23ll2^ |. 
but not in the ordered phase. Although models of this 
type do not account the effects of orbital and spin de- 
grees of freedom, important to the physics of real metals, 
it captures the basic physics of a quantum phase transi- 
tion in a metallic Fermi system. As we shall see below, it 
describes the phenomenon of 'local quantum criticality' 
observed in many experiments. 

The bosonization of dense Fermi systems, relativistic 
or not, in d > 1 space dimensions requires the intro- 
duction of local coordinate frames near the FS to de- 
scribe the fluctuations of this smooth curved quantum 
object. Much as in the conventional Landau theory of 
the Fermi liquid^, the fluctuations of the FS describe 
the quantum diffusion of particle-hole wave packets on 
the FS 17] . We therefore introduce a set of N patches 
of width A = 2irkp/N ~ \Jkp\ given by the curvature 
of the surface and cutoff with N — > oo in the renormal- 
ization group (RG) limit of A — > 0. Thus, all correc- 
tions of order 1/N are uninteresting within this frame- 
work provided they correspond to irrelevant operators. 
A key feature of the bosonization treatment is that the 
corrections to the linear dispersion (locally normal to the 
FS) are described by non-linear operators (in the bosons) 
which are irrelevant in the Landau phase and at the fixed 
point theory. Processes which involve large(er) momen- 
tum transfers tangent to the FS are described by both 
intra-patch and inter-patch interactions. Contrary to 
some claims I^Et)]. bosonization treats the effects of the 
curvature of the FS exactly. 

Introducing local coordinates for patch S, located at 
Fermi wavevector ks, we find that this phase space anal- 
ysis allows us to expand the fermion operator, ip(x), in 
terms of a set of fermion species derived from each patch: 
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(3) 

In terms of the local Fermi velocity vs at patch S, a 
vector locally normal to the FS, the dispersion is 



Zq = vs ■ q + P(v s ■ q} 2 + j(vs ■ q) c 



(4) 



where — ctc^ is the fermion occupation number oper- 
ator, N(0)f 2 (q) — F 2 /{1 + n\F 2 \q 2 ) is an interaction of 
range y/n, 
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is the Fourier transform of the nematic order parameter 
field, the quadrupolar density, and £ Q is the quasiparti- 
cle energy which is kept up to cubic order in \q\ — Uf- 
This model was proposed in Ref. p| to study a transition 
to a nematic Fermi fluid which in mean field occurs for 



where \vg\ — vf and \q\ -C kp- Due to this dispersion, 
our fluctuating FS can be thought of as N x oo fermion 
species each with its own conserved charge due to the 
U(l) Nxo ° symmetry of the Landau theory and of Eq.Q: 



ips(x»,x t ) -> e ies{ - Xt) i} S (x n ,x t ) 



(5) 



where (x n ,xt) are local orthogonal patch coordinates 
with x n parallel to vs- Though we have presented this 
within the language of our local coordinates, this charge 
conservation for each point of the FS is a well known 
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property of the Landau Hamiltonian, e.g. Eq.Q, quite 
independent of the choice of coordinates on the FS. 

A basic consequence of Eq.JSJ) is a simplification of the 
fermion correlation functions. If one creates a fermion of 
type (S, x' t = 0) at (x' n = 0,t' = 0), one must destroy 
another fermion of the same type at (x n , t) to have any 
overlap with the ground state. Hence, only one species 
may contribute to the single particle fermion Green func- 
tion. After defining a windowing function W&(xt) that 
vanishes for |xt|A > 1 and attains the value of unity at 
Xf = 0, the correlator becomes 

G f (x,t) = 

£ ^^e tt **» (G\T T ^ s (x n , 0, r)4(0, 0, 0)|G) (6) 

s 

where x n = k$ ■ x and x t _L x n . Following standard 
bosonization techniques we compute G f (x,t) for T > 0: 



G f (x,t) =22 Tr G F{s) (x n ,T)Z s (x n ,T) 

s 

(7) 

where G° F ^(x n , r) is the 'patch' free Green function and 



In Z s (x n ,T) = 



y 1 - 

Z-^ R 



d 2 q 



(2tt) 2 ^ (3 



} i(k n x n —uj n T) 



V s R £ A (q,iUn) 
(iuj n -vs ■ qf 



(8) 



with V<?g A the effective RPA interaction. 

The limit N — > oo can be obtained. Consider the case 
when Zs is independent of the patch index S (true for a 
circular FS) and decays or oscillates slower than e lkpXn 
as x n ^ oo . Converting the sum over patches to an 
integral, and using steepest descents, we find that since 
the integral is peaked around 6 = cos -1 (/cs • x/k F r) = 0, 
we get the angular average ( e ik F rcos6 } — J (k F r), and 



G F {r,t) ~ G%(r,t)Z 8=0 (r,t) 



(9) 



with G° F (r,t) ~ cos(kpr)/r 3 / 2 ; the weighting function 
Wn(xt) and sum over patches have dropped out. 

Near the nematic QCP of Q, Vg'g A is directly related 
to the propagator of the order parameter field. Although 
for angular momentum channels greater than s-wave, the 
order parameter theory contains a variety of collective 
modes, it has been shown that InZs is always domi- 
nated by the softest of these modes, a z = 3 over damped 
mode very similar to that occurring in Hertz's original 
treatment of the ferromagnetic problem. 

Following Refs.P, 0, 0], and using standard results 
of the well known nematic-isotropic classical phase 
transition}^, we construct the qualitative phase dia- 
gram for the Nematic quantum phase transition of a 



Fermi fluid. There are clearly three regimes of interest. 
Two of them are disordered and described by 



yBPA 

v s,s 



{q, i^n 



1/N(0) 



(10) 



nq z 



where A = A(S,T) with 5=1 + 1/F 2 , the distance to 
the mean field position of nematic QCP at F% = — 1. In 
the FL regime, A is dominated by the scale 8 and we 
obtain Aw 5 — 5c where 5c arises from finite corrections 
to the gaussian theory due to irrelevant interactions. In 
the quantum critical regime, A is dominated by temper- 
ature and we must resort to a careful renormalization 
group argument to understand its behavior Q- We find 
that the temperature dependence of A follows the law 
A = — 7TTn(T), where 7 is the coupling to the danger- 
ous irrelevant quartic interactions of the order parameter 
theory, which in the model of Eq.Q is the cubic term 
in the dispersion relation, 

HEil Eq.Q. In the nematic 
phase, the Goldstone mode dictates that A = 0, 



yRPA 
V S,S 



1/N(0) 



(11) 



qv F 



r 28 s 



In contrast to Ea. (|10[l . it depends on the patch index S. 

We now compute G f (x,t) as a function of T and <5, 
near the QCP. We are interested in the behavior of the 
amplitude Zs, discussed at T = in Ref.[l^|. For a FL, 



const 
(12) 



const 



Z s (x n ,t = 0)U„^oo = Z <l, Z s (0,i)|t^oo 

At the QCP, the leading behavior is 

Z s (xn,t = OJU^oo ~ e-( a "/ R ) 1/3 , Z S (0,t)\t^ o 

(13) 

There is a similar non Fermi liquid behavior in the ne- 
matic phase, except that there are four special points at 
sin 29s = where the leading behavior is that of a FL. 

For T > 0, the uj n = contribution to Zs in Eq.JHJ) 
dominates for T K <C T -C T F , T K = v F /kBK, and 



InZ, 



I static 



N(0)k 



1 



d 2 q cos(q n x n 
2ir 2 v 2 F ql(q 2 + A/k) 



(14) 



k B T 



N(0) 



^<log(C/x„) 0<x„«e 



X n > ^ 



(15) 



where ^ = y/ k/A. On the FL side, ^ is quite small and 
we obtain exponentially decaying spacial correlations. In 
the quantum critical regime, £ is very large, set only by 
the irrelevant interactions and temperature. Finally, in 
the nematic phase £ = 00 for T > 0, and Zs — * 0, kept 
finite only in a finite sized system (except once again at 
the four special points). At equal positions (x n = 0), 
the Lo n = contribution vanishes, and the equal position 
auto-correlation function is well behaved in the infrared. 
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Near the QCP, In Zs (x n , 0) obeys the scaling law: 

lnZ s (x n , 0, S, T) = 6- 1 / 3 In Z s (te n , 0, «5/6 2/3 , T/6) (16) 

where all length scales are measured in units of the 
interaction range \fk. Upon the change of variables 
q' n = q n /b, q' t = qt/b 1 / 3 , as b -> 0, 9 -> q'6 1 / 3 , 

(17) 

For T = <5 = we find, 
In Z s (x„, 0,0,0) = |x„| 1 /3 lnZs( sign (a; n ), 0,0,0) (18) 

by letting b = |a; n | _1 . Since In Zs(±l, 0, 0, 0) < 0, we 
find Zs(x n , 0, 0, 0) ~ e a. | a= | 1 / 3 ag ex p ec t ec i_ From 

Z = lim Z s (x„,0, <5,T) = lim {n kp - q - n kF+q ) (19) 

and Eq.JHJ, we find Z = at the QCP. As 5 -> at 
T = 0, by letting b = 5 3/2 we get 

Z = exp (\nZ s (oo,0,l,0)/y/s\ (20) 

Thus, since lnZs(oo, 0, 1, 0) < 0, Z vanishes with an 
essential singularity as 6 — * 0. Thus, this scaling form, 
Z$(x n , 0) clearly exhibits the breakdown of the FL as the 
QCP is approached. 

Let us complete our discussion of scaling by going to 
T > 0. Letting A = £~ 2 and b = T in ifTHfl . we obtain: 

hxZ s (x n ,0,r 2 ,T) =T" 1 / 3 In Z s (s„T, 0, T 2 ^ 2/3 , l) 

(21) 

Since we only need Zg for T > 0, we may evaluate it 
again using the uj n = mode (now in the scaling limit): 

\nZ s {x n ,{),C 2 ,T) = -A\x n T\i/k F (22) 

so that holding x n T fixed, this diverges like 
(— TlnT) -1 / 2 as T — > and therefore it is not a 
scaling function of x n T. Since the free fermion Green 
function in J7J asymptotically behaves as e~ n ' XnT \ we 
see that the above form dominates over the free fermion 
value for ^_ > it, or for T < T g = T K / (gk F n), where we 
have reintroduced n with T K — vf/ y^ks- Hence, for 
T 3> T g , the free fermion result dominates and scaling 
occurs in x n T. For T -C T g , scaling in x n T breaks down. 

In this letter, we demonstrated non-perturbatively the 
phenomenon of local quantum criticality near the metal- 
lic charge nematic quantum phase transition. It ap- 
pears to be a robust behavior of many quantum crit- 
ical systems, as, in addition to the nematic quantum 
critical point, it also applies to the critical behavior of 
quantum dimer models [l2j, metalic antiferromagnets and 
ferromagnets[bH . While the order parameter correla- 
tions, £ are large in the quantum critical regime, kept 



finite only by the dangerous irrelevant quartic interac- 
tions, the fermionic correlations in this regime are cor- 
respondingly short ranged, kept non-zero by the same 
dangerous irrelevant quartic interactions. Indeed, in the 
nematic phase, which has an infinite correlation length, 
we have shown that the corresponding fermionic correla- 
tions vanish in the thermodynamic limit except at four 
special points on the FS. An extension of this analysis 
to other quantum phase transitions in metallic systems 
results in the same basic conclusion, that the 'local quan- 
tum criticality' of the fermions is independent of z and 
it is due to the oj n = static contribution to the corre- 
lation function. We found that while the single particle 
density of states, N{oj) remains finite throughout the en- 
tire phase diagram, the qp residue Z (the jump in the 
fermion occupation number) vanishes for T <T g both in 
the quantum critical regime and in the nematic phase, 
which suggests that a FS may not be seen in a de Haas- 
van Alphcn experiment. 
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